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We investigate the transport properties for a Luttinger liquid coupled to two identical Bose-Einstein conden- 
sation reservoirs. Using the approach of equation of motion for the Green function of the system, we find that 
the distance between the two resonant transmission probability peaks of the system is determined by the bosonic 
interaction strengths, and the sharpness of these resonant peaks is mainly determined by the Rabi frequency and 
phase of the Bose-Einstein condensation reservoir. These results for the proposed system involving a Luttinger 
liquid may build a bridge between the controling transport properties of cold atom in atom physics and the 
interacting boson transport in low-dimensional condensed matter physics. 
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I. INTRODUCTION 

The physics of ultracold one-dimensional (ID) Bose system 
is very different from that of ordinary three-dimensional cold 
gases JH01- The possibility of Bose-Einstein condensation 
(BEC) in one dimension has been discussed for the noninter- 
acting Bose gas yfl. The interaction between bosons plays an 
essential role due to the strong constraint in phase space in ID 
case. Monien et al. 0] have shown that a trapped quasi-one- 
dimensional system of interacting bosons under the experi- 
mental conditions can be described by a Luttinger liquid (LL) 
Hamiltonian. As is known that the low-energy physics of ID 
single channel conductors can not be described by the Fermi 
liquid theory if the particle-particle interactions are taken into 
account |5|]. Such system falls into the so-called LL regime. 
However, unlike the Fermi liquids, the LL liquids may also 
include ID interacting bosonic systems. Bosonic systems can 
display fermion-like properties and vice versa (a, |8|] . One 
well-known example in the field of cold atoms is the behavior 
of the Tonks-Girardeau gas [9], where the bosons interact so 
strongly that they effectively behave as free fermions. 

In LL model theory, the main assumption is the lineariza- 
tion of the free-particle dispersion relation near the eigenen- 
ergy points of the system. Fermionic systems which are be- 
lieved to be described by the LL model include quasi-one- 
dimensional organic metals IToll . quantum wires fulfill , and 
edge states in the quantum Hall system [13]. The actual sys- 
tem considered is finite length and are attached to two identi- 
cal reservoirs at its end points. This is to say, these systems 
are always embedded in a three-dimensional matrix. So they 
will show a crossover to three-dimensional behavior at low 
temperature, while the trapped ID Bose gas would provide a 
clean testing ground for the concept of LL model. 

In nanoelectronics the control of electron quantum wires or 
quantum dots is performed by the biased conducting leads at- 
tached to them. But in nanobosonics the role of the "leads" 
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is replaced by the finite superfluid reservoirs (given particle 
numbers) which can be coupled to a particular atom by opti- 
cal transitions. With regard to this field, the dynamics of an 
atomic quantum dot coupled to a BEC reservoir via laser tran- 
sitions has also been studied recently lfl4l[T5ll . 

In this paper, we consider a system consisting of a LL cou- 
pled to two identical BEC reservoirs. The bosons in the BEC 
reservoirs are confined in a shallow trap, while the atom in the 
LL is confined in a very tight potential. Atoms in both the LL 
and BEC reservoirs correspond to the different internal atomic 
states connected by Raman transition with Rabi frequency Q. 
and detuning 5. Using the approach of standard equation of 
motion for Green function (GF), we investigate the frequency- 
dependent transport properties for this system. Our results 
show that the distance between the two resonant transmission 
probability peaks is determined by the interaction strengths, 
and the sharpness of the resonant peak is mainly determined 
by the Rabi frequency and the phase of the BEC reservoir. The 
results for the proposed system involving a LL may build a 
bridge between the atomic transport in atomic physics and the 
interacting electron transport in low-dimensional condensed 
matter physics. 

II. MODEL 

The total Hamiltonian of the system consists of three parts, 
i.e., 

H = J] H a +H ce „+H T , (1) 

o=L,R 

where H a is the Hamiltonian for the isolated left or right BEC 
reservoir, H cen is the Hamiltonian of the isolated LL, and Hj 
is the Hamiltonian describing the transfer of a particle from 
the BEC to the LL. 

The starting point for the calculations of H a (a = L,R) is 
the Hamiltonian 

H a = f rfr[-^(r)ivV»(r) + V(r)^(r)^„(r)] 
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r )Mr )if, a (r), (2) 



where if/dr) and ifr a (r) are annihilation and creation operators 
for bosons in the BEC reservoir respectively, m is the atomic 
mass, V(r) is the potential confining bosons system in a trap, 
and U(r - r ) is the interaction potential between two parti- 
cles in the BEC reservoir. (We have adopted the unit of h-\ 
throughout this paper.) To take into account the quantum fluc- 
tuations of the state in which all the atoms are condensed in a 
single quantum state, the operator iff a (r) can be represented in 
the form ^i a {r)- -\/No<f>o+6i//a(r), where No is the particle num- 
ber in the zero-momentum state, O is the wave function of the 
condensed state, and 6iff a (r) denotes the fluctuation operator 
of momentum A: + 0, i.e. the excitation above the ground state. 
Within the Bogoliubov approach one assumes that 6\j/ a (j) is 
small and retains in the interaction all terms which have two 
powers of ifr a (r) or ty a (r). This is equivalent to including terms 
which are no more than quadratic in 8ift a (r) or 6tf>l(r). 
Performing the Fourier transformation 



6ifr a (r) = 2^ Off. 
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replacing 6ifr a (r) and 6if/ a (r) by a a ^ and a , for k + 0, we 
obtain 



H a = ^[(<e° Jc + Ui(k))(al k a ayk + al_ k a a - k ) 

+ U 2 (k)(al k a ] ak + OajcOa-k)], 



(4) 



where e°, is the single particle energy, and U\{k) (U2(k)) is 
the Fourier transformation of U(r - r ). Here the operators 
a atk and a ak are destroy and create bosons in the state with 
momentum k satisfy the usual Bose commutation relations. 

We take the effect when two atoms are close to each 
other into account by using the effective interaction, and the 
Hartree-Fock terms are both equal to «of/o in which «o is 
the number density of the BEC, Uq the contact interaction in 
Hatree-Fock approximation, so the Hamiltonian for the iso- 
lated left or right BEC reservoir reads lfl6ll 



H a = ^ e «,k + noU Q )(al k a^ k + a^_ k a a -k) 
+ n Ua(al k al_ k + a a , k a a - k )], 



(5) 



where the single particle energy e° , = k 2 /(2m), and a^ ak (a a , k ) 
is the creation (annihilation) operator of the bosons in the left 
or right reservoir. Note that the prime on the sum indicates 
that it is to be taken only over one half of momentum space. 

Considering the Raman detuning between LL and the reser- 
voir and the phase of the reservoir, the Hamiltonian of the 
isolated left or right BEC reservoir, i.e., the first term H a 
(a = L,R) of Eq. (1), is given by 



Ha = X ^ 6a ' k + \^ a i,k a ^ k + a l,-k a " -*) 



+ (Aa I ,a\ , + A*a a , k a a - k )], 



where e a , k = e , + 6 a in which 5 a is the Raman detuning be- 
tween LL and the reservoir a; A - noUne^" in which <p a is the 
phase of the reservoir a. In the absence of the driver (8l=Sr), 
the particles can also pass through the system because of the 
phase difference between two BEC reservoirs. 

For the Bose gas in a cylindrical symmetric trap confined 
to the z axis by a tight trapping potential in x-y plane, if the 
extension L of the trap in the z direction is much larger than its 
radius R and the temperature is much lower than the energy of 
the lowest radial excitation, the ground state is described by a 
LL 0]. The starting point for the calculations of H cen is also 
the Hamiltonian 

+ \f dr f dr>J(r)^(r')l/(r-rVi(rVi(r)j(7) 

where i//\(r) and if/[(r) are annihilation and creation operators 
for bosons in the LL respectively, V (r) is the potential con- 
fining bosons system in a trap, and U (r — r ) is the interaction 
potential between two particles in the LL. Through the same 
procedures as above, the operator i//i(r) can be represented in 

the form \fr\{r)= .JN^(f>' +Stf/i(r). And the fluctuation operator 



6ifn(r) = Yjb q 



(8) 



where the operators b q and bl that destroy and create bosons 
in the state with momentum q satisfy the usual Bose commu- 
tation relations 

= V. [&„2y]=0, [Aj,^] = 0. (9) 

Using the Bogoliubov approach and replacing 6tfri(r) and 
8tff^(r) by b q and b q for q + 0, we obtain 

Hcen = Yj [(e ° + Ul(qm l b " + b -9 b -" ) 
1(1*0) 



+ U 2 {q){blb\ + b q b- q )}, 



(10) 



where e° is the single particle energy, and U\(q) {Uiiq)) is the 
Fourier transformation of U (r — r ). 

It is known that in the LL model, the main assumptions are: 
(1) the linearization of the dispersion relation; (2) only small 
momenta exchanges included. In terms of the two assump- 
tions, the Hamiltonian of the isolated LL, i.e., the second term 
H ce „ of Eq. (1), is given by |T3l 

Hcen = ^[{vdq + kO + ^Kbpfq + bl^b-g) 



q>0 

V 2 + t 



(11) 



(6) 



Here the single particle energy e q -VL(q + k£) because of the 
linearization of the dispersion relation. In Hamiltonian (11), 
bt (b q ) is the creation (annihilation) operator of the bosons in 
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the LL, vi is the eigen-velocity in the channel, ki is the eigen- 
wavevector, and V) (j = 2, 4) is the interaction potential when 
q ~ in which V 2 represents the interaction potential between 
the left- and right-moving boson branches, while V4 repre- 
sents the interaction potential within a momentum branch. 

Note that the LL arisen in our system describes Bose sys- 
tem, so the operators b q and b\ correspond to the destruction 
or creation of an individual particle (i.e., boson). However, 
when the LL describes the interacting electrons in one dimen- 
sion, an individual particle is a fermion. By means of the 
bosonization technique, we can also write the Hamiltonian in 
terms of boson operators. But the operators b q and b' q are lin- 
ear combination of the density fluctuations p q = £ ctck+q, so 

k 

they conserve the number of fermion particle and do not cor- 
respond to the destruction or creation of an individual particle. 

The Bose field can also be describe by its density-phase 
representation: ipB{r)- ^p{r)e' ei - r) . Expanding in small fluctu- 
ations of the phase 69 and the density dp around the saddle- 
point solution, i//B(r)= Vpo + 6p(r)e i[9(,+S9(r>] . The density fluc- 
tuation operator dp and the phase fluctuation operator 59 form 
a pair of conjugate operators [69(z), 5p(z )]=iS(z — Z ). With 
the same approximation as the equations of motion in Ref [4], 
we can also express the Hamiltonian of the isolated LL as (Eq. 
(10) in Ref. [4]) 



H a 



s 



dz 



J-(d z 59) + —dp 1 
Lm Lp L 



(12) 



where p is the number of particles per unit length, m is the 
atomic mass, and k is the compressibility. 

The starting point for the calculations of Hj which de- 
scribes the transfer of a particle from the BEC to the LL, is 
the Hamiltonian 



h t = n L ^(nWi(n) + o^IWjfte) 



(13) 



where Q, a (a = L, K) is Rabi frequency. Here we have assumed 
that the atom in the LL is coupled to atoms in the reservoir a 
via Raman transition with Rabi frequency. Using the Bogoli- 
ubov approach and replacing 5iff a (r) and 6ip' a (r) by a a ^ and 
a , for A: ^ 0, we obtain 

a,k 

H T = nL^e^X/iMn) 

k 

+ Q* <5-A"[(r 2 ) J] a R , k e ikn + h.c, (14) 



And replacing 6ifri(r) and 5i[r\(r) by b q and b\ for q + 0, the 
Hamiltonian (14) becomes 



H T = ^\Q a a[ k b q + Q.* a b\a a , k ). 

k,a,q 



(16) 



Because the operator a} . correspond to a creation of a particle 
in the BEC reservoir and the operator b q correspond to a de- 
struction of a particle in the LL, a terms such as b\a a ^ would 
thus correspond to a destruction of a particle in the BEC and 
a creation of particles in the LL. 

However, for a fermionic LL, if there is transfer of a particle 
from the BEC reservoirs to the fermionic LL, it must have a 
different form, 

H T = Q L i/r[(r 1 )i/r F (r 1 ) + Q*^(r 2 )(AR(r 2 ) 

+ n^(n¥i(n) + ( 17 ) 

where i[fF=T* c q e' q r is the Fermi annihilation operator. Using 

1 

the bosonization technique, the Fermi annihilation operator 
can be written as [18] 

~ J] <! fe*M + *M ~ exp(iYje iq "'b q + h.c.)), (18) 

p=±\ q 

where §(r) and <p(r) obey the commutation relations 
[<p{r\ 0(r')]=[#(r), #(/)]=0 and [<p(r), §(r')]=isgn(r - r')/2. 

Note that the total Hamiltonian H in our system is equal to 
the sum of Eq. (6), Eq. (11) and Eq. (16). 



III. FORMULATION 

In terms of the Heisenberg equation of motion, the current 
of the reservoir a can be written as 

dN a k 

Ia,k(t) = (-£*■) = -Wafi.m) 

at 

= -i(Yj&aaljp q )+h.c.), (19) 

kq 

where A^,* = £ a t k a a,k is the total number operator for the 

k 

boson in the reservoir a. Defining a 2x2 GF G < ia (f), then the 
current of the reservoir a becomes 



I a ,k(t) = -/^(CV3^ ff (f))ii +h.c, 



(20) 



kq 



where h.c. denotes the Hermitian conjugate. In terms of the 



operator i//i(r)= ^N' (f> + differ) = y/p + 6p(r, fj e m+ee ^ ] , if 
replacing 5ij/\{r) by e' se (r), the Hamiltonian (14) becomes 



Ht 



-ik-n t jse{n) 

u L,k 



Z 



ik'i^ 



+ h.c. 



(15) 



where (G < ka (t)) 1 1 = (a a k b q ) is the element in the first row and 
the first column of GF G < ka (t). 
Similarly, 

= -i Yj^ G< qM(t))22 + h.c, (21) 

kq 

where (G^ ka (t))22 is the element in the second row and the 
second column of GF G < , (f)- Since the current is conserved, 

q,ka - 7 ' 
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the currents of the bosons with momentum k and —k are equal, 

I.e., I a k — Ia,-k- 

Using the theorem of analytic continuation, we have 



(22) 



where G < ' r is 2x2 lessor/retarded GF of LL with cou- 
pling between the LL and the reservoir, while g <! k is 2x2 
lessor/advanced GF of the isolated BEC reservoir, respec- 
tively. The GF gi JE) in Fourier space is given by 



8:jE) = K/E) ~ 8*, k (E)lfa(E), 



(23) 



where f a (E) (a = L,R) is the Bose distribution function. And 
based on the Landauer-Buttiker formula lflS^l20ll2~ill . the cur- 
rent in the Fourier space can be written as 



r°° dE 



[f L (E)-f R (E)]T(E\ 



(24) 



where E is the energy of the incident boson from the reservoir 
a. And with the help of Keldysh equation: G < = G 1 "L < G a , the 
transmission probability T(E) is solved as 

T(E) = 4Tr[r L (E-6 L ) YG r {E)L r R {E-d R ) Vg»(£)]. (25) 



1 
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Here E < ^ r is lessor/retarded self energy, respectively. The re- 
tarded self energy is given by V = Y_ "E' a where U a is retarded 

self energy of reservoir a and defined as 



(26) 



where 



Q„ 
-Q.1 



(27) 



and g r , is the retarded GF of the isolated left or right BEC 
reservoir and defined as 



{{a a ,k(t), a ak (t))) {{a a ,k(t), a a - k (t)}) 
{{al_ k (t), al k (i))) <{o^(0. a a ,-kit)}) / 



(28) 



In terms of the equation of motion for the GF, g r . in the 
Fourier space can be written as 



= 1 / E + e ak + |A| -A 

'■ k E 2 -e 2 ak -2\A\ eak \ -A* -E + £ „k + \A\ 



(29) 



Defining g a =^ g r a t> by straightforward calculation, we ob- 

k 

tain 



Pa 



E -A 
-A* -E 



r d{e aA + |A|) 

J (£ + ; +)2 - (^ + |AD 2 + 



/^t + lA| o \ r 



|A| 2 

d{6a,k + |A|) 



•*0+) 2 -(e^ + |A|) 2 +|A| 2 
(30) 

Here we have changed Y into an integral f c/ec, kP a with the 

it J 
help of the density of the states in the BEC reservoir p a . The 

second term of Eq. (30) vanishes, because in the second term 

bothg^jj and g' all are odd functions of e^£+|A|. So we finally 

obtain 



, _ I E -A \ f d{e a , k H 



+_|A|) 

+ |A|) 2 + |A| 2 



■(31) 



In the following, calculating the intergral by using the residual 
theorem, Eq. (31) can be reduced as 



—ivn 



y/E 2 + |A| 2 



2 -A* -E 



(32) 



with v = 1 for E > and v = — 1 otherwise. 

Inserting Eqs. (27) and (32) into (26), we obtain 



ivT a 



1 



E A 

2 V^ 2 + IAI 2 I A * ~ E '* 



(33) 



where the linewidth function Y a =2n\£l 2 a \p a . Under the so 
called wide-band approximation, the self energy of the lead 
is not sensitive to the energy and can be taken as a constant 
independent of the energy E. The non-diagonal term in the 
expression of the self energy A=| A|e'^ is | A| when the phase of 
the BEC reservoir 0=0. 

In Eq. (25), G r (E) and G"(E) denote the Fourier transforms 
of the GF G r (t) and G"(f) respectively. G r can be obtained by 
Dyson equation in matrix form 



G r - G r Q + G r rG r , 



(34) 



where GJ, is the retarded GF of the isolated LL which is de- 
fined as 

Gq(M) = -iO(t-t) 



({b q (t),bl(t')}) ({b q (i),b-Jt)}) 



({bl g (t),b q (t)}) ({bL q (t),b- q (t)}) 
In terms of the equation of motion for the GF, we can obtain 



(35) 
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G r = 



1 



1 



1 



E L - (v L q{- + l)/2 + v L k L Y + (v L q(- - \)I2) 1 



E + v Lq a + l)/2 + v L k L - v ^(i-l)/2 



-vmAjt - l)/2 



-£ + v L? (4 + l)/2 + v L £ L J 



(36) 



where the parameter g is the strength of the interaction which 
is defined as g — (1 + V/(7rvz,))~ 1/2 . Here we have assumed that 
V2-V4-V. This definition follows that of the fermions. The 
LL parameters g also can be extracted from the Lieb-Liniger 
equation 11221 12311 . For repulsive bosons, g-\ corresponds to 
the hard-core limit, while g>\ for repulsion, with g— >oo in 
the limit of weak interactions. In the case of fermion, non- 
interacting fermion corresponds to g-1 and repulsive interac- 
tion corresponds to g < I. 



IV. QUANTUM TRANSPORT 

In the following we show some numerical examples calcu- 
lated according to Eq. (25) for the transport properties of this 
system with the experimental parameters: 02411 for 87 Rb, T-l 
nK and the eigenenergy of the LL £^=2.0 kHz. By analyzing 
the form of G r Q for the isolated LL, because q/k^Q or qjkjj^l, 
and the energy is equal the sum of the excitation energy and 
the eigenenergy of the LL, there should appear peaks near the 
eigenenergy and near three times of the eigenenergy in the 
transmission probability versus the energy of the incident bo- 
son. And peaks near E/Ei-l or near E/El-3 will evidently 
differ from the resonant peak in the case of quantum dot cou- 
pled to BEC reservoirs. Because q has a range of values, this 
will open some new channels for transmission. Here we will 
consider the symmetric case, i. e., F£=r«=r/2. 

Fig. 1 illustrates the transmission probability T as a func- 
tion of the incident boson energy E (in units of Ei) where the 
parameters are 5l-5r-Q and 0=0. For the system with fixed 
|A|=27r x 0.41 kHz and g=10, Fig. 1(a) shows the dependence 
of two different Rabi frequency £2=0.02 kHz (red solid line) 
and 0.03 kHz (blue dashed line) on the transmission probabil- 
ity, respectively. From this figure we can find that the reso- 
nant peak becomes wider as the Rabi frequency increases. It 
is because |Q| 2 oc F, and F describes how well the reservoir 
is in contact with the LL. The larger linewidth function corre- 
sponds to the stronger coupled case, and the stronger coupling 
corresponds to the wider the resonant peak. Fig. 1(b) shows 
the result of the transmission probability T versus E with fixed 
Rabi frequency Q=0.02 kHz and |A|=2ttx0.41 kHz for two 
different interaction strengths, where the red solid line for g-2 
and the blue dashed line for g=10, respectively. The distance 
between the two resonant peaks become smaller when the 
particle-particle interaction parameter g is larger. From this 
figure, we can conclude that the interaction parameter plays 
an important role on the relative position of the two resonant 
peaks. Fig. 1(c) illustrates T as a function of E with fixed 
Rabi frequency 0=0.02 kHz and fixed interaction strength 
g=10 for two different |A|=0.27r kHz (red solid line) and 0.82^ 




E 



FIG. 1 : (Color Online) The transmission probability T as a function 
of incident boson energy E (in units of E L , E L =2.0 kHz) where the 
parameters are 6l=Sr=0 and 0=0, (a) for two different Rabi frequen- 
cies £2=0.02 kHz (red solid line) and 0.03 kHz (blue dashed line) 
with |A| = 2^x0.41 kHz and g=l0, (b) for two different interaction 
strengths g=2 (red solid line) and 10 (blue dashed line) with £2=0.02 
kHz and |A| = 2^x0.41 kHz, and (c) for two different |A|=0.2tt kHz 
(red solid line) and 0.827T kHz (blue dashed line) with £2=0.02 kHz 
and #=10, respectively. 



kHz (blue dashed line), respectively. From this figure, we can 
not find more visible difference between the two cases, which 
show that the width and height of the peaks are not sensitive 
to |A|. 

Fig. 2 shows the dependence of three different phases of 
BEC reservoir on the transmission probability with fixed Rabi 
frequency £2=0.02 kHz, |A[=2^x0.41 kHz and g=10, where 
the blue line for 0=0 or (p=2n, the red dashed line for 0=7r/2 
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E 



FIG. 2: (Color Online) The transmission probability T as a function 
of the incident boson energy E (in units of E L , E L -2.0 kHz), wtih the 
parameters of S L =S R =0, n=0.02 kHz, g=10 and |A| = 2^x0.41 kHz, 
the blue solid line for <p=Q or 2n, the red dashed line for (j> = n/2 or 
3n/2, and the green dashed line for (f> = n, respectively. 



or (f>=3n/2 and the green dashed line for (p-n, respectively. 
From this figure we can see that in the range of to n the 
resonant peak becomes sharper as the phase increases, while 
in the range of n to 2n the resonant peak becomes wider as 
the phase increases. Fig. 2 is very similar to Fig. 1(a), which 
makes clear that, through varying the off-diagonal term of the 
self energy, the phase of the BEC reservoir play the similar 
role as the Rabi frequency on the transmission probability T. 

There may be possible experimental realizations for our 
system. Firstly, two BEC reservoirs can be realized in current 
experiments with atomic gases. Secondly, Luttinger liquid 
arisen in our systems can be realized in current experiments 
with atomic gases. With the current technology there ap- 
pears no difficulty in making transverse frequency « ± >100w,, 
where a> z is the longitudinal frequency. In such limit, one can 



produce atomic gases with all the atoms lying in the lowest 
harmonic oscillator state in the x-y plane, leaving the motion 
along z (the only degree of freedom). The system then be- 
haves like a ID Bose gas. For steeper magnetic traps, a>± ~ 
50 kHz, particle densities of p ~ 10 4 particle/cm, and assum- 
ing a scattering length of 1 10 ag for Rb, it should be possible 
to observe the LL behavior |4J]. 



V. CONCLUSION 

In conclusion, using the equation of motion for Green func- 
tion, we have investigated the transport properties for a Lut- 
tinger liquid coupled to two identical Bose-Einstein condensa- 
tion reservoirs. It is demonstrated how the transmission prob- 
ability is determined by Rabi frequency, interaction strength, 
|A|, and phase of the BEC reservoir, respectively. We have 
found that the distance between the two resonant transmission 
probability peaks is determined by the interaction strengths, 
while the sharpness of the resonant peak is mainly determined 
by the Rabi frequency and phase of the reservoir. The fur- 
ther theoretical investigation on taking into account impurity, 
spin or other interactions are worthy to be carry out. These 
results for the proposed system involving a LL may be useful 
to control transport properties of cold atom. 
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